Neutron specific heat in the crust of neutron stars from the nuclear band theory 
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The inner crust of neutron stars, formed of a crystal lattice of nuclear clusters immersed in a sea 
of unbound neutrons, may be the unique example of periodic nuclear systems. We have calculated 
the neutron specific heat in the shallow part of the crust using the band theory of solids with Skyrme 
nucleon-nucleon interactions. We have also tested the validity of various approximations. We have 
found that the neutron specific heat is well described by that of a Fermi gas, while the motion of 
the unbound neutrons is strongly affected by the nuclear lattice. These apparently contradictory 
results are explained by the particular properties of the neutron Fermi surface. 

PACS numbers: 

Keywords: neutron star crust, neutron star cooling, band theory, specific heat, dynamical effective mass 



Neutron stars are born in the gravitational core col- 
lapse of massive stars. With about one or twice the mass 
of the Sun compressed inside a radius of only 10 kilo- 
meters or so, neutron stars are among the most compact 
objects in the Universe The outer layers of the star, at 
densities below the neutron drip threshold pnd — 4x 10 11 
g.cm~ 3 , are formed of a solid Coulomb lattice of neutron 
rich nuclei coexisting with a degenerate gas of relativis- 
tic electrons. The inner crust, at density above pnd and 
below the crust-core transition density, which is at about 
half the density p = 2.8 x 10 14 g.cm -3 inside atomic nu- 
clei, is permeated by a sea of unbound neutrons, which 
could be superfiuid in some layers 0. 

Many observed neutron star phenomena are intimately 
related to the physics of the crust. For instance, X-ray 
bursts in Low-Mass X-ray Binaries are associated with 
thermonuclear burning in the accreting neutron star en- 
velope 0| . In the subgroup of soft X-ray transients, accre- 
tion outbursts are followed by long period of quiescence 
during which the accretion rate is essentially zero. In 
some cases, the period of accretion can last long enough 
for the crust to be heated out of equilibrium with the 
core. The thermal relaxation during the quiescent state 
has been recently monitored for KS 1731—260 and for 
MXB 1659—29 after an accretion episode of 12.5 and 2.5 
years respectively Q. These observations thus provide 
information on the thermal properties of neutron star 
crusts [H, 0] • The thermal relaxation of the crust could 
also be potentially observed in very young isolated neu- 
tron stars, 10 — 100 years after their formation 

The diffusion of heat in the inner crust is mainly gov- 
erned by the thermal properties of the shallow layers, ow- 
ing to their very low thermal diffusivity [9j. One of the 
key parameters is the neutron specific heat. In cooling 
simulations of neutron stars, the neutron specific heat is 
generally approximated by that of uniform neutron mat- 
ter. At low enough temperatures, neutrons are predicted 
to become superfiuid but the exact value of the critical 



temperature T c still remains uncertain (2j. Moreover, the 
presence of nuclear clusters is likely to modify the critical 
temperature. For instance, in Ref. :4], it has been shown 
that the pairing field in the dilute neutron gas completely 
vanishes just after the drip point. 

The effects of superfluidity on the neutron specific heat 
arc usually incorporated in the following way [Io| : 



paired 



(T)=TZ(T/T c )c v (T) 



(1) 



where the renormalisation factor 1Z(T/T C ) accounts for 
the effects of superfluidity and cv(T) stands for the spe- 
cific heat calculated in the absence of superfluidity. In 
the above expression, the presence of the nuclear clus- 
ters is ignored. Their effects have been studied by solv- 
ing the Hartree-Fock-Bogoliubov equations for the nu- 
cleons [ll], [ijj and by including the contribution of col- 
lective modes These calculations have been car- 
ried out within the spherical Wigner-Seitz approxima- 
tion. However the validity of this approximation has 
been recently discussed Pointing out the analogy 
between unbound neutrons in neutron star crusts and 
conduction electrons in ordinary metals, the band theory 
of solids have been applied and adapted to the nuclear 
scale [lj| [l6| . It has shown for the first time the large 
effects of Bragg scattering on the motion of the unbound 
neutrons. One may expect that thermodynamical prop- 
erties like the specific heat could also be modified. 

In this paper, we investigate the effects of the solid 
crust on the neutron specific heat, within the framework 
of the nuclear band theory (Sec. [TJ) . Given the present 
uncertainties in the pairing problem, even in the simpler 
case of pure neutron matter Q , we will focus on the nor- 
mal part cy(T) only. Results in the shallow layers are 
discussed in Sec. |TT] and are compared to several approx- 
imations used in the literature. 
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I. A MICROSCOPIC MODEL OF THE 
NEUTRON STAR CRUST 

Following the standard assumptions, we consider a 
body centered cubic crystal with only one nuclear cluster 
per lattice site [l[ . We have taken the composition (N Z) 
of the clusters calculated by Negele and Vautherin [17| • 
But the nucleon distributions have been recalculated at 
each temperature T by solving the finite temperature 
Hartree-Fock equations with the Skyrme SLy4 effective 
nucleon- nucleon interaction [l8|. Typical temperatures 
of interest for cooling isolated neutron stars and for ac- 
creting neutron stars lie in the range between 10 to few 
10 9 K @. 

Fully self-consistent calculations within the band the- 
ory are computationally very expensive. We have thus 
solved the Hartree-Fock equations in two steps [3]. 
First, we have determined self-consistently the nucleon 
distributions and mean fields in the Wigner-Seitz ap- 
proximation with mixed Neumann-Dirichlct boundary 
conditions. We have removed the spurious fluctuations 
of the neutron density by averaging in the interstitial 
region [17]. In a second stage, we have solved the 
Schrodingcr equation with Bloch boundary conditions 
using these self-consistent mean fields. We have fixed 
the lattice spacing so that the volume of the Wigner- 
Seitz polyhedron is equal to the volume of the spherical 
Wigner-Seitz cell of radius i? C eii- We have employed the 
Linearized Augmented Plane Wave method described in 
details in Ref. [lj|. As shown in Ref. 14], the one-body 
spin-orbit potential has a negligible effect on the neutron 
level ordering and was therefore neglected. 

The specific heat (per unit volume) is defined as 



Cy(T) 



dU 
df 



T 



ds_ 

df 



(2) 



(taking the Boltzmann's constant k-Q = 1) where U is 
the total internal energy density, S the entropy density 
and the partial derivatives are evaluated at constant vol- 
ume V . We have found that the latter expression is nu- 
merically more convenient. The entropy density of the 



TABLE I: For the three shallow layers considered in this pa- 
per, are given: the total density (p), the number of protons 
and neutrons (Z, N respectively), the nuclear mass number 
A (at T = 0), the radius of the cell (iiceii) and the neutron 
gas density at T=0 (p™). 



p [g.cm -3 ] 


Z 


N 


A 


R ccU [fm] 


P° [fm- 


- 3 ] 


6.69 x 10 11 


40 


160 


133 


49.24 


1.3 10" 


-4 


1.00 x 10 12 


40 


210 


141 


46.33 


2.6 10" 


-4 


1.47 x 10 12 


40 


280 


140 


44.30 


4.9 10" 
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unbound neutrons is given by 
d 3 k 



(2tt 3 ) 



fak In fak + (1 - fak) hl(l - f a k) 



(3) 



where a is the band index, k is the Bloch wave vector and 
f a k is the Fermi-Dirac distribution. The latter depends 
on the temperature and the neutron chemical potential 
Hn which is determined from the total density p„ of un- 
bound neutrons by 



E 



d 3 fc 

(27r3) /afc - 



(4) 



Integrations have been carried out using the special point 
method (see [Hj]). Note that the number of unbound neu- 
trons per cluster is determined by the shape of the mean 
fields, which vary with temperature. Consequently p 1 ^ 
can also depend on the temperature. In particular, with 
increasing temperature the most loosely bound neutrons 
may be excited into the delocalized states. 



II. RESULTS AND DISCUSSIONS 

We have considered three different layers in the shallow 
region of the inner crust. Their properties are summa- 
rized in Table [U The neutron chemical potential /i n (T) 
and the entropy density S(T), have been calculated for 
a set of temperatures between T = 10 keV and T = 100 
keV. We have then evaluated the specific heat accord- 
ing to Eq. ([2]), by interpolating the entropy density S(T) 
with cubic splines and numerically differentiating. 

For the range of temperatures considered here, T <C 
£f, where £f is the Fermi energy. The specific heat is 
then approximately given by (ref. [lj|, p47) 



cv(T)~— g(e F )T, 



(5) 



where <?(£f) is the density of states at the Fermi level. 
This expression is valid provided i) the density of un- 
bound neutrons p„ and £f are independent of the tem- 
perature, ii) the density of states g(e) is sufficiently 
smooth. Pertaining to the first point, we have found 
numerically that p% varies very little so that the neu- 
tron evaporation/condensation phenomenon mcntionncd 
in Sect. [Jean be ignored. In a previous paper, we have 
shown that the density of states in the shallow layers of 
the crust is far from being smooth at a scale of few keV 
(see Fig. 6 of Ref. [lj])- However, since we consider ther- 
mal energies of order 10 keV or more, the variations of 
the density of states with the energy are expected to be 
much smaller at this scale as noticed in Ref. [l4| • As can 
be seen in Fig. [TJ the neutron specific heat varies almost 
linearly with T as expected from Eq. ([5]). 

In order to understand more qualitatively our results, 
we have performed a comparison to several approxima- 
tions that have been used in the literature. 
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A. Validity of various approximations 

The simplest approximation to the neutron specific 
heat is to neglect the presence of the nuclear clusters. 
The specific heat of a Fermi gas is given by 



c v /T [fm^MeV 1 ] 



FG 



2/3 m®T 



,V3 



(6) 



where m® is the Skyrme effective mass evaluated at T = 
for the density of the neutron gas p„. Note that for 
the layers we considered, to® ~ m n . As shown in Fig. [TJ 
the specific heat given by © is very close to the exact 
result from the band theory. 

This striking result can be understood using 
the semi-classical Extended- Thomas-Fermi expansion 
method [2(J. At lowest order in fi 2 , the neutron specific 
heat is given by 



TF 



(T) 



7T\2/3 T 



d 3 r 



™>K(r) 1/3 . (7) 



We have found that higher order corrections in U are 
negligible. The specific heat ([7]) shown in Fig. [1] for the 
three selected layers, is very close to the results of band 
theory. The small effects of the nuclear lattice can now be 
easily explained by the small volume fraction occupied by 
the clusters (typically ~ 10~ 2 — 10 ~ 3 ). Indeed in the limit 
of uniform neutron matter, the Thomas-Fermi specific 
heat ([7|) reduces to Eq. ((6|). A more refined explanation 
within the full quantum mechanical framework will be 
given in Sect. Ill 51 

We have also computed the specific heat with the sin- 
gle particle energies e a obtained in the Wigner-Seitz ap- 
proximation. Varying the temperature, we have found a 
small redistribution of the neutron energies e a close to 
the Fermi energy. Due to the very large number of neu- 
tron states, a small redistribution of energies can have 
a large impact on the total energy and on the entropy. 
However, it has been found numerically that the entropy 
being a smooth function of the occupied and unoccupied 
states, is much less affected by this redistribution than 
the total energy. We thus have calculated the specific 
heat using 



ds ws 



dT 



The entropy density is given by 



5' 



ws 



cell 



/«ln/ a + (l-/ a )ln(l-/ ) 



(8) 



(9) 
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FIG. 1: (color online) Neutron specific heat divided by the 
temperature in the three shallow layers calculated with the 
band theory (solid thick line) and with different approxima- 
tions: free Fermi gas c^ G ([6| (dotted line), Extended- Thomas- 
Fermi Cy F © (dashed line), Wigner-Seitz Cy S © (solid line 
with filled squares) and de Gennes c^ G (|10p (solid line with 
filled triangle). 



where g a is the degeneracy of the state a. The specific 
heat ([H]) is shown in Fig. [1] Despite non-linear fluctu- 
ations due to the redistribution of neutron states, the 
results are fairly close to those of the band theory. The 
best agreement is obtained for the layer with the lowest 
density, as expected from the domain of validity of the 
Wigner-Seitz method [l4l |. 



We have checked the validity of another approxima- 
tion, suggested by de Gennes [2l[ and applied in Ref. [IH , 
consisting in differentiating the entropy ([9]) with respect 
to the temperature but assuming that the energies e a and 
the chemical potential /j, n (as well as the pairing gaps in 
superfluid systems) are independent of T. The specific 
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heat is then given by 



(10) 



As can be seen in Fig. [TJ this approximation is valid at 
very low temperature where the energies and the chemi- 
cal potential could be taken independent of the tempera- 
ture but it becomes less and less reliable with increasing 
temperature. 

From the comparison of various approximations to the 
neutron specific heat, we have found that the Fermi gas 
model reproduces very well the results of the band the- 
ory, while more sophisticated approaches like the Wigner- 
Seitz or the de Gennes approximations are less accu- 
rate. Does it mean that the unbound neutrons are really 
"free"? 



B. Are unbound neutrons really "free"? 

The present results on the specific heat together with 
Eq. ([5]) indicate that the average density of states are 
close to that of a Fermi gas (as already noticed in a pre- 
vious paper, see Fig. 6 and 9 of Ref. This conclusion 
is at first sight surprising since the mean field potential is 
very deep inside the clusters (about -70 MeV) . On general 
grounds, one expects the effects of the nuclear lattice to 
be negligible whenever the Fermi wavelength of the un- 
bound neutrons Af = 27r/fcp with fcp = {iit 2 p^) 1 ^ 3, , is 
much larger than the lattice spacing. However this con- 
dition is not satisfied for the layers we considered. Using 
the densities p„ given in Table Q] we find Af — 40.09, 
31.82 and 25.76 fm for R ceU = 49.24, 46.33 and 44.3 fm, 
respectively. One therefore expects the unbound neu- 
trons to be strongly scattered by the clusters. 

If the unbound neutrons were really free as suggested 
by the results on the specific heat, this would imply that 
their motion is unaffected by the lattice. It is well-known 
in solid state physics that the interactions between the 
conduction electrons and the ionic lattice can be taken 
into account through renormalizating the mass of the 
electrons. The motion of an electron in a solid, with a 
wave vector k in a band a, can thus be characterized with 
a dynamical effective mass tensor defined by (ref. [l9l |. 
p228) 



1 

m* 



1 



-(e) 
'ak 



h 2 dkidkj 



(11) 



where fci, i = 1,2,3 are the Cartesian components of 
k and the electron energy. This concept has been 
extended to low-energy thermal neutrons propagating in 
crystals and the corresponding effective mass has been 
experimentally measured in silicon [22| . In the context of 
neutron star crust, since the number of unbound neutrons 
per lattice site can be very large, it is more appropriate 



to introduce an average effective mass to* defined by [1 

d 3 A; 1 d 2 e a k 



K 



It, 



F (2tt) 3 h 2 dhdh 



(12) 



where the integral is taken over all occupied states. This 
dynamical effective mass to* can be equivalently ex- 
pressed as 



1 



1 



d 2 u 



-Y 

£ ^r 1 dpnidpr, 



(13) 



where U is the energy density of the moving neutrons 
in the crust frame and p n is the average neutron mo- 
mentum. This effective mass has implications for neu- 
tron star dynamics. For instance it has been shown 
that a large enough effective mass can trigger a Kelvin- 
Helmholtz instability which might explain the origin of 
pulsar glitches [231 ]. 

The dynamical effective mass to* has been calcu- 
lated for the considered layers, using the same numerical 
method as in Ref. [16[ . Results are shown in Table [Til 
If the unbound neutrons were free we would have found 
to*/to„ = 1. However, the dynamical effective mass is 
much larger than the bare neutron mass to*/to„ > 1 in- 
dicating that the interactions between the neutrons and 
the clusters are very strong. In the following, we will 
explain the apparent contradiction with the results ob- 
tained previously for the specific heat. 

At the beginning of Sect. [TT1 we have shown that the 
specific heat varies almost linearly with the temperature 
(see Fig. [IJ, as expected from Eq. (J3J). The density of 
states #(ef) appearing in Eq. (j3J) can be expressed as an 
integral over the Fermi surface <Sf, 



1 



E 



dS 

s F Kit I 



(14) 



where v a k — h~~ Vfc£ Q fc is the group velocity of the un- 
bound neutrons. The presence of nuclear clusters leads 
to the existence of resonances. As a result, the energy 
bands of unbound states may be locally flat in k-space 
thus distorting the Fermi surface The deformation 
can be estimated by the ratio £f of the area of the de- 
formed surface to that of the unperturbed Fermi sphere. 
Therefore, we have iSp = £F47rfcp, where by definition 
£f = 1 in the absence of clusters. Numerical calcula- 
tions show that the nuclear lattice reduces the area of 
the Fermi surface by more than a factor of 2 in the shal- 
lowest layer (see Table [TTJ . However since the specific 
heat is driven by that of a Fermi gas, we can infer from 
Eq. (TT4")) that in average, the group velocity is changed 
by the same factor £p compared to that of the Fermi gas. 

Likewise the coefficient K, (fT2"|) can also be written as 
an integral over the Fermi surface 



1 



3(2tt) 3 ^ 



Yf Kk\dS. (15) 
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From the previous discussion, it follows that this coeffi- 
cient is therefore approximately changed by a factor £p. 
Eq. (fl~2f then implies that m* w rn„/£ F . Inspecting Ta- 
ble |T] shows that the predicted relation between m* and 
£f is only roughly satisfied because v a k varies on the 
Fermi surface. Nevertheless, this analysis in terms of the 
topology of the Fermi surface explains why the presence 
of the nuclear clusters has such a strong impact on the 
motion of the unbound neutrons but not on their specific 
heat. 



The results on the specific heat might suggest that the 
unbound neutrons are not affected by the presence of the 
nuclear clusters. It is however not true since the dynam- 
ical effective mass of the same unbound neutrons is very 
different from the bare one unvealing strong interactions 
with the periodic lattice. This apparent paradox can be 
explained by the fact that the neutron Fermi surface area 
and the neutron group velocity are both reduced by the 
same factor in the presence of clusters, leading to a nearly 
unchanged density of states and specific heat. 



III. CONCLUSIONS 

Modelling the recently observed thermal relaxation of 
neutron star crusts requires the knowledge of their ther- 
mal properties. In this paper, we have computed the 
specific heat of the neutron ocean permeating the inner 
crust at densities below 1.5 x 10 12 g.cm" 3 , by applying 
the band theory of solids with the Skyrme SLy4 nucleon- 
nucleon interaction. We have compared the results ob- 
tained using different approximations. We have found 
that for temperatures T — 10 7 — 10 9 K relevant to neu- 
tron stars, the neutron specific heat is essentially given 
by that of a Fermi gas. The Thomas-Fermi expression ([7]) 
yields nearly undistinguishable results. The specific heat 
calculated in the Wigner-Seitz approximation agrees rea- 
sonably well, while the de Gennes approximation ([TO]) 
leads to specific heats a factor of 2-3 smaller than those 
obtained in the band theory. 



TABLE II: For the three shallow layers considered in this 
paper, are given: the total density (p), the dynamical effective 
mass (m*/m n ) and the reduction factor of the Fermi surface 
area (£f). 



P [g-cm -3 ] 


m* n /m n 






6.69 x 10 11 


4.0 


0.44 


0.77 


1.00 x 10 12 


3.6 


0.49 


0.86 


1.47 x 10 12 


3.2 


0.52 


0.87 
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